Coherent edge mixing and interferometry in quantum Hall bilayers 



Stefano Roddaro, 1 ' 2 Luca Chirolli, 3 Fabio Taddei, 1 Marco Polini, 1 and Vittorio Giovannetti 3 ^ 

1 NEST, Istituto Nanoscienze-CNR and Scuola Normale Superiore, 1-56126 Pisa, Italy 
2 Istituto Officina dei Materiali CNR, Laboratorio TASC, Basovizza, 1-34149 Trieste, Italy 
3 NEST, Scuola Normale Superiore and Istituto Nanoscienze-CNR, 1-56126 Pisa, Italy 

We discuss the implementation of a beam splitter for electron waves in a quantum Hall bilayer. 
Our architecture exploits inter-layer tunneling to mix edge states belonging to different layers. We 
discuss the basic working principle of the proposed coherent edge mixer, possible interferometric 
implementations based on existing semiconductor-heterojunction technologies, and advantages with 
respect to canonical quantum Hall interferometers based on quantum point contacts. 
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I. INTRODUCTION 

Chiral edge states living at the bound ary of a two- 
dimensional (2D) quantum Hall (QH) phased constitute 
a fascinating playground both for the investigation of fun- 
damental properties of one-dimensional electron liquids 3 
and the implementation of innovative devices. 

A number of novel QH electron interferometers have 
been demonstrated in recent years, opening a new win- 
dow of investigation on coherent electron transport in 
solid-state devices. In particular, edge "beams" have 
been exploited for the realization of a variety of elec- 
tronic interferometers reproducing the Mach - Zchn- 
dei^Hll, Fabry - Perot 7 and Hanbury Brown - TwissP^l 
schemes adopted in optics. These can have an important 
impact both on the fundamental investigation of quan- 
tum transport phenomena of electrons in solids and as 
possible implementations for quantum computing^. In 
these circuits, mixing between edge states has so far been 
achieved using beam splitters (BSs) based on quantum 
point contacts (QPCs). Fascinating but still puzzling 
phenomena have been highlighted in these devices^, in 
particular in relation to finite-bias visibility and edge re- 
construction phenomena. While the QPC approach has 
proven successful, the intrinsic geometry of this BS imple- 
mentation makes it necessary to adopt non-simply con- 
nected 2D electron gases (EGs) and limits the complexity 
and size of the achievable 2D QH circuits. In addition, 
such BS structures represent a potentially comp lex c ir- 
cuit element displaying non-linear characteristics^^ 3 -' as 
well as, in some cases, fractional substructures^ which 
can have an impact on the overall behavior of the 2D 
QH circuit. Alternative interesting device schemes exist 
and are based on tunneling between co-propagating edge 
modes^ED. 

Here, we discuss a different paradigm to edge-beam 
interferometry, which is based on the exploitation of a 
QH bilayeA This is a system composed of two closely- 
spaced 2DEGs. We assume that each of the two 2DEGs 
is in the QH regime and has edge states. Our BS for 
electron waves is based on inter-layer tunneling (Fig. [I]), 
i.e. hybridization between the two 2DEGs. We show 
that momentum-conserving tunneling between the two 




FIG. f . (Color online) A top view of a quantum Hall bilayer 
interferometer (QHBI). In this cartoon, the QHBI is obtained 
by employing two overlapping orthogonal two-dimensional 
electron gases ( "top" and "bottom" conducting strips of width 
W) in the integer quantum Hall regime. The two subsystems 
are vertically separated by a distance h. We assume that only 
inter-layer tunneling couples the two subsystems. In this con- 
figuration, transport takes place only via the quantum Hall 
edge states of the two subsystems. Inter-layer tunneling is 
effectively active only at the edge crossing points, labeled by 
A, B, C, and D. 



2DEGs can be used to scatter an edge mode localized 
in one layer into an edge mode localized into the other 
in a controllable way. In particular, when edge modes 
in the two 2DEGs cross at a finite angle, the scattering 
process effectively involves on ly a regio n of the order of 
the magnetic length, Ib = \Jhc/{eB), around the geo- 
metrical crossing point of the edge modes. Using this 
approach, many interferometric schemes can be achieved 
by employing relatively simple gating geometries. 

This Article work is organized as follows. In Section |ll| 
we describe a basic interferometric scheme, which can be 
obtained by overlapping two orthogonal two-dimensional 



2 



electron gases in a bilayer system. Section III describes 
the formal scattering problem ruling the interferometer 
behavior. Section [TV] discusses the adopted scattering 
matrix for the beam splitter and analyzes the details of 
its microscopic origin. Finally, in Sect. |Vl| we summarize 
our main findings and draw our conclusions. 



II. INTERFEROMETER LAYOUT 

We consider a quantum Hall bilayer interferometer 
(QHBI) constituted by two two-dimensional electron 
gases (2DEGs) (the bottom layer, B, and the top layer, 
T) extending in the x, y plane and separated by a dis- 
tance h along the vertical direction z. Both 2DEGs are 
assumed to be in the integer quantum Hall regime, in- 
duced by the presence of a (uniform) quantizing mag- 
netic field B pointing along the z-axis, i.e. B = Bz. 
We assume that the two subsystems are coupled only 
by a uniform tunneling term of strength Asas- While 
Asas is mostly determined by the heterostructure de- 
sign (height of the barrier, inter-layer separation h, etc) it 
can also be tuned to some extent by an in-plane magnetic 
fielcP^I. In principle, the two 2DEGs are also coupled by 
electron-electron interactions. In this Article, however, 
we neglect these effects and discuss the working principle 
of our QHBI at the single-particle level. The inclusion 
of inter-layer electron-electron interactions is definitely 
challenging and expected to be responsible for an inter- 
esting phenomenology, which is, however, well beyond 
the scope of the present Article. 

We further assume that the electrostatic potentials 
that are responsible for the lateral confinement within 
each layer, are "layer-dependent" and characterized by 
mutually orthogonal "longitudinal" directions. Specif- 
ically, introducing the in-plane coordinate vector r = 
xx + yy, we take the confining potential Vb(t) in the 
B layer to be translationally invariant with respect to x, 
while the one Vx(r) in the T layer to be translationally 
invariant with respect to y, i.e. 



V B (r) = V B (y) 



Vr(r) = Vi(x) 



(1) 



The specific functional form of Ve(y) and Vr(x) will be 
fixed later. Under this condition and at low energies, 
charge transport is dominated by single-particle chiral- 
edge modes propagating along the x axis in the B layer, 
and along the y axis in the T layer, granting the setup 
the form of two rectangular Hall bars of width W, which 
are crossing perpendicularly as schematized in Fig. [T] 
Analogous configurations have been realized experimen- 
tally in semiconductor-heterojunction double quantum 
wells thanks to top and bottom gating and to inter- 
layer screening effects^"2iJ. These allow depleting one of 
the two quantum wells without substantially altering the 
other and thus can lead to independent carrier-density 
profiles in T and B layers. Significant tunneling can be 
obtained for suitably designed barriers. 



A. Model Hamiltonian 

The Hamiltonian of our system consists in the sum 
of two terms, H = Ho + Htun, which describe respec- 
tively the free evolution of the electrons in the two lay- 
ers and the inter-layer tunneling coupling. Introducing 
Fermionic field operators ^ CT (r) that satisfy canonical 
anti-commutations rules^, the first contribution can be 
expressed as 



Ho = J2 I d2r &Ar)H a (r)i> a (r) 

er— T,B 



(2) 



where 



2m* 



-m r + -A{r) +V a {r) . (3) 
c 



Here m* is the electron band mass (m* ~ 0.067 m e , 
for example, for GaAs, where m e is the electron mass 
in vacuum) and A(r) is the vector potential associated 
with the uniform magnetic field B = Bz. We work in 
the symmetric gauge, A(r) = (Bx r)/2. The inter-layer 
tunneling term can be written a£l 



H t 



*SAS 



*t( r )*B(r) +H.c. 



(4) 



where Asas is the symmetric-to-antisymmetric tunneling 
gap. Values of Asas ranging from <C 1 fic\ -'' up to 
the meV scale^ have been demonstrated in literature by 
tuning the barrier design and the in-plane magnetic field. 

Equations ^ and Q can be casted in a more compact 
form by expressing them in terms of the eigenstates of the 
"unperturbed" problem (i.e. Asas = 0). To this end, 
we introduce the eigenvalues e a ,n.k of the single-particle 
Hamiltonians H a (r) and the corresponding eigenfunc- 
tions <t>„, n ,kX r ) : 

0B,«,*(r) = e-™y/^^= XB.nAy) , 

V Z7T 

„iky 

<h,n,k(r) = e™y^^= X T,n,k{x) , (5) 

where k € K is the eigenvalue of the magnetic transla- 
tion operator along the x (y) axis for the B (T) layer, 
while n > is a discrete index which labels the Landau 
levels. The functions Xcr,n.k{ua) determine the transverse 
structure of the propagating modes ([5| and are the eigen- 
functions of the transverse Hamiltonian: 



h 2 d 2 

2m* dv?„ 



\u a - M CTjfc ) 2 + V„{u a ) 



Here w c = eB/(m*c) is the cyclotron frequency, Ub,, 
and u„ is the "transverse" coordinate: 



(6) 



~UT,k 



"B ■ 



ub = y 



(7) 



Since the functions 4> a ,n,k{ r ) form a complete orthonor- 
mal seP^, we can use them to expand the operators 
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*& CT (r), obtaining the identities 

OO p 

*<x( r ) = ^2 dk <^,r> : fc( r ) C<r,n,, 

rt J 



n=0 ' 



Ccr.n.k 



d'r 0* fe (r) * ff (r) 



(8) 



(9) 



where c aini k are the Fermionic annihilation operators as- 
sociated with the chiral-edges modes of the unperturbed 
Hamiltonian "Ho that satisfy the eigenmode equation 



Analogously, we can write 



OO p 
td t> <n — n " 



^B.Tn^O 



and 



(10) 



(ii) 



Q2) 



with the tunnel matrix elements defined as 
Asas 



f 



n' ,k' :n,k 



A 



SAS c -jkk 

An 



d2r 4>T,n>,k'( r )<l>li,n,k( r ) 

•F n',n(k , fc) , (13) 



where T n ',n is a form factor that describes the overlap be- 
tween transverse wavefunctions residing in different lay- 
ers: 



■Tv.nOfc',*:) = / dxdy e -*(*+*'<l)(v-"l)/3> 

X XT,n'/ S '(a ; )XB,n,fe(y) • (14) 



III. THE EIGENVALUE PROBLEM IN THE 
PRESENCE OF TUNNELING 

To characterize how inter-layer tunneling affects the 
transport properties of the system at hand we need to 
study the eigenmode equation 



-e7(e) 



(15) 



which plays the role of Eq. ( 10 1 when tunneling is taken 



into account. More explicitly, Eq. (15) can be written 



for the components of the pseudospinor wavefunction 
V (e3 (r) = (^ £) (r),4 e) (r)) T , i.e. 



H B (r)^\r) 



^(r) = e4 6) (r) 



(16) 



which make it explicit that the eigenfunctions of the full 
Hamiltonian are delocalized in the two layers due to the 



tunneling coupling Asas- Both ( 15 ) and ( 16 1 can also be 
casted as 

(e - eB,n,fc)41*: = J2j dk ' *n',fe';n,fc «T,n',fc' > 
n' 

TV 

(17) 

where a a ,n,k{^) are complex coefficients relating 7(e) to 
its unperturbed counterpart c a ,n,k and the components 
of the spinor ip^ e '(r) to the functions (f> a -,n,k,{ r ) defined 
in Eq. ([5]), i.e. 

ct— T,B n ^ 

^ £) (r-) = E / *41 fe ^,W- (19) 

n J 

At this point it is clear that the total field operator 
^( r ) = J2<7^v( r ) can l* e expanded in the basis of the 
eigenvectors 7(e) as 



*(r) = ^ ^ (e) (r) 7 (e) , 



(20) 



where e compactly labels all quantum numbers of the 



problem ( 16 ) 



The eigenmode equation ( |16| defines a scattering pro- 
cess in which an edge beam described by the unperturbed 
energy eigenstatc ([5]) and propagating in the B layer, say, 
is partially transmitted into the T layer. The associated 
probability amplitudes strongly depend on the specific 
choice of the confining potentials ([I) and are, in general, 
rather difficult to compute due to their complex func- 
tional dependence on e CT! „,fe, t n ',k>;n,k, and 4> a .n,k{r). 

The analysis greatly simplifies in the weak tunneling 
regime where Asas can be treated as a small pertur- 
bation compared to all other energy scales in the prob- 
lem by employing the Born approximation. Consider, 
for instance, the event corresponding to elastic scatter- 
ing (eB.n.k = eT,n',fe') from the unperturbed energy eigen- 
state 4>B,n,k(f) (which describes an electron in the n-th 
Landau level, propagating with momentum hk along the 
x-axis on the B layer) to the state 4>t ,n' ,k' i r ) ■ Elemen- 
tary manipulations of the ordinary equations of scatter- 
ing theorySS in the Born approximation yield the follow- 
ing expression for the associated probability amplitude: 



S[<j>B,n,k -> <fa,n',k'] — (-27Ti) 

d 2 r 



Asas f j2 ( r ) <pB,n,k( r ) 



where 



zAsas e lkk ' l B 

2h y/VT,n',k'VB,n,k 

1 de an k 



Va.n,k 



h dk 



F n , in (k',k), (21) 



(22) 
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is the group velocity of the mode 4> a .n,k{ r ) an( i where 
we have used Eq. ( 13 1 to express the result in terms of 
the form factor ( 14 ) . Due to the perturbative nature 
of Eq. ( pT| ), its validity is restricted only to those cases 
where the modulus of S[<pB.n,k ~^ <t>T,n',k'] is small (we 
shall provide momentarily a more precise statement on 
this). Going beyond this regime is typically extremely 
challenging. 

Surprisingly, though, our scattering problem admits an 
explicit analytical solution in the special — but typically 
valid — scenario of smooth confinement. In this limit the 
change of the confining potentials ([I]) over a magnetic 
length is assumed to be negligible with respect to the 
cyclotron gap, i.e. 



dV a (u a ) 



di 



(23) 



This assumption well describes the physics of edge states 
defined by electrostatic gates and has two main conse- 
quences, both extremely useful in simplifying the scat- 
tering problem. Specifically, 



i) it implies a tunneling Hamiltonian (12) where 
modes c CT , n! fc, c a ' . n ' ,w with n ^ n' are effectively de- 



coupled; smooth confining potentials obeying ( 23 ) 



do not mix edge beams living in the two layers and 
corresponding to different Landau-level indices. 

ii) under proper conditions, it allows the linearization 
of the dispersion relation of the unperturbed energy 
eigenvalues e CTn fe. 



These properties originate from the fact that Eq. (23) 



permits to well approximate the transverse wavefunctions 
X<y,n,k entering Eqs. (pi) with the bulk Landau level eigen- 
functions^, i.e. with (properly translated) eigenstates of 
a harmonic oscillator of frequency lo c and zero-point fluc- 
tuation £b, 

-V2, 



XB,n,k(y) ' 
XT,n : k( X ) 



1/2 II 



U\)ll B ) , 
U B )/i B ) , 



(24) 



where tp n (£) 



-« 2 / 2 ff n (£), H n (0 being the 



th Hermite polynomial. At the level of the unperturbed 
energy eigenvalues of the system this implies that one can 
write them as a harmonic contribution plus a correction 
associated with the confinement potential, i.e. 



e^n.fc ~ huj c {n + 1/2) + V a (kl%) 



(25) 



which can then be turned into a linear expression in k as 
required by property ii) by properly limiting the inter- 
val of k which enters the problem [see the following for 
details]. Property i) instead follows by observing that 
within the approximation ( 24 ) the form factor ( 14 ) of 



the system becomes diagonal in n and independent from 
the momenta k and k', i.e. 



2tt£ 



B f 



'2tt £ 



B %" 



(26) 



For future reference, notice that, in writing the second 
identity, the operator algebra of the harmonic oscilla- 
tor^ has been adopted to express J" n ',„ in terms of the 
system Fock states \n) and \n'). Eq. (26) implies that 



the Landau-level index n retains its validity as approxi- 
mate quantum number even in the presence of tunneling. 
Hence we can drop the summation over n in Eq. ( 19 ) 
looking for eigenmodes of the form 
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(e) = [ dk ( a l%r^n, k , (27) 

a=T,B J 

i>i%(r)= f dk ai% ^ n , k Xr) , (28) 



where now ct^ k solve the following simplified system of 
coupled equations: 



(e) _ M)"A S 



2V2TT 



I Jkk'tL (e) 



j7 t ikk £ R 

ak e B a 



T,n,k' ' 



(29) 



/ \ ( e ) i n A 

(e T ,n,fc - eja T ,n,fe 



2v / 2tt 



dk' e***» a^ k , 



The smooth confinement condition ( 23 ) turns out to 



be useful also to better clarify the limit of validity of the 
weak tunneling expression (21). Indeed, let us focus on 



the special case in which the two edge states involved in 
the scattering process belong to the lowest Landau levels 
of their respective layers {i.e. n= nf = 0) and propagate 
with the same momentum hk under smooth confinement 
conditions. According to |24] ), in this case the transverse 
wavefunctions Xc.o.k can be approximated by Gaussian 
amplitude distributions with variance Iq. Their overlap 
vanishes in an area of radius ~ £b around the geometric 
crossing point of the associated classical skipping orbits, 
yielding a form factor proportional to £b — see Eq. p6| — 
and a corresponding a scattering amplitude (21) with a 
modulus that scales as Asas^b/(^f)- Here, vp is the 
associated group velocity of the two modes which we as- 
sume to be identical. The condition for the validity of 
the perturbative approach can thus be casted in terms of 
the following inequality: 



VF A S as 



(30) 



This admits a simple physical interpretation in terms of 
the ratio between the time £b/vf spent by an electron 
crossing the active tunneling region of size £b at a speed 
Vf, and the time S/AgAS that is necessary to tunnel from 
one layer to the other one. 



A. A specific example: linear confinement 

In this Section we illustrate more explicitly the notion 
of "smooth confinement" by discussing a specific exam- 
ple, which is amenable to a fully-analytical treatment: 
the case of linear confinementPsl. 
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FIG. 2. (Color online) Interlayer edge mixing in the limit of 
a linear confinemen t po tential, which obeys the smooth con- 
finement condition ( 23 1 . Edge modes propagate west to east 



in the bottom layer and south to north in the top layer. Inter- 
layer tunneling is only active in a region of size w Ib around 
the geometric crossing point of the edge modes. The math- 
ematical analysis of the resulting scattering model leads to 
an ideal beam-splitter, closely resembling a semi-transparent 
mirror, with a transparency that depends monotonically on 
the strength of inter-layer tunneling Asas ■ The dependence of 
the confining potentials V a {u a ) on the transverse coordinate 
u a is illustrated in the right panels. 



Let us assume that the potentials ([I]) have a linear 
dependence upon the transverse coordinate u a in each 
layer: 



V B (y) = eEy 



Vr(x) 



-eEx 



(31) 



E being the intensity of an applied uniform electric field. 

With this choice, the unperturbed eigenenergies of the 
system acquire a linear dispersion on momentum hk, 



1 



Hvpk — -m*Vp 



£<j,n,k = huj c (n + 1/2) 
characterized by a constant group velocity 

E 



v F 



eEl% 



(32) 



(33) 



Note that vp coincides with the classical expression for 
the drift velocity in crossed uniform electric and magnetic 
field. 

Furthermore, the transverse wavefunctions Xan,k 
can be easily written in terms of properly translated 
harmonic-oscillator eigenstates: 

XB,n,k(v) = ^B 1/2( Pn((y - k£g)/£ B + Co) , 

XT, n ,k{x) = r B 1,2 y n {{x + kl\)lt B - Co) , (34) 
where we have introduced the dimensionless parameter 

eEt B 



6 



hlOr 



(35) 



Such functions describe chiral edges modes which, for 
k > 0, propagate from west to east in the B layer, and 



from south to north in the T layer (see Fig. [2|. The 
corresponding form factor ( 14 1 is independent of k and 
k' and equal to 



n' ,r 



= V2ir 



£ B i n ' e< /2 (n'\D (ie iff/4 £ ) \n) . (36) 



Here, as in Eq. (26 1, we have used the operator algebra 
of the harmonic oscillator^! to find the result expressed 
by the last equality. The displacement operator D(a) is 
defined b^ 



D(a) = exp (aa* — a* a) 



(37) 



where a (a') is the harmonic-oscillator destruction (cre- 
ation) operator. 

The parameter Co gauges the smoothness of the linear 
potential with respect to the cyclotron gap. The cri- 
terion Co *C 1 is indeed the smooth-confinement condi- 



tion for the case of the linear-confinement model (31 1. 
(Notice in particular that in the limit Co ~> Eq. (34) 



yields Eq. (24). Equation (36) makes it explicit that for 
Co — > 0, the coupling between different Landau levels 
decreases exponentially with their distance, with a lead- 



ing term which is proportional to Co 



\n — n\ 



In particular, 



Eq. (36) reduces to the diagonal expression (26) in the 
limit Co 0, implying that tunneling only couples edge 
beams with the same n and allowing us to look for eigen- 
mode solutions of the form (28). Exploiting this fact and 



Eq. (32), the corresponding eigenmode equation (29) can 



be finally casted in the following form: 



[K F k\ a B n k - 



Jdk' 



u T,n,k' 



(38) 



[4" } - k] a, 



(e) 
T,n 



7 % 
V2k 



J dk' e 



-ikk'l 



B a 



(e) 

B,n,k' ' 



where we have introduced the quantities, 
Asas ^b 



k 



2 hv-p 

(„) _ e — hw c (n + 1/2) + m*v F /2 



(39) 
(40) 



As discussed in the Appendix [A] Eq. (38) admits ana- 
lytical solutions of the form 



*B,n.fc 



= e ikk^h% f(°\ k i B -k F n> i B )e B 



x T,n.k 



(41) 

k ( ; i h B )) i B , 

(42) 



with 



-7 [c 0(/c) + d 0(-«)] \np /«, (43) 
-7 [a 0(-«) + b 0(k)] \k\-^ 2 /k. (44) 



6 



Here 0(k) is the Heaviside step function, while a, 6, c and 
d are complex coefficients fixed by the boundary condi- 
tions. These impose the following linear relationships 
between c, d and a, b: 



7 



d = 



2tt 

7 



rH7 2 ) 
r(-z 7 2 ) 



a e 



-7T772 



(45) 



where T(- • • ) is the Euler Gamma function. By replac- 
ing (41) and (|42| into (28 1 we finally get the eigen- 



functions. For the sake of simplicity, we here report 
only those associated with the lowest Landau level (i.e. 
n = 0). Specifically, introducing the dimensionless pa- 



rameter £p = kp £b and the variables Sx = x 
Sy = y — £f^b, and z = (Sy + i5x)/£b, one has 



V'b.oM = -7 



B.O.fci 



>(r) e 



x [c A+(-«-f ) -d A+(z + £o)] 



V't.oO) = ~7 ^o^lW 



(46) 
(47) 



where 
and 



x [-a A_ + f ) + & A_ (-iz - f )] , 
• o fe (0) ( r ) are * ne e ig eniu nctions (151 at Asas = 



A ± (z) 



r(± i7 2 )e z2 / 4 ZWz) , 



D K (z) being the parabolic cylinder special function. In- 
terestingly, mathematically similar solutions have been 
independently obtained in a very recent worlpU analyz- 
ing the QH effect in a 2DEG subject to a potential oc xy. 

For \z\ ^> 1 they admit the following asymptotic ex- 
pansion 



A±(z)~r(± i7 2 



(49) 



which can be used to study the asymptotic behavior of 
the solutions (46) and ( [47] ) . In particular, from this it 
follows that in the limited x/Ib | ^> I, the B-layer com- 
ponent behaves as 



'2-n 



At) e l? p e -n 2 in|5x/te| 



x [a Q(Sx) + be(-Sx)] 



(50) 



where the constraint ( 45 ) was employed in simplifying 



the expression. Similarly, for \Syj£s\ > 1 we have 
x [cQ(Sy)+de(-6y)} . 



(51) 



Equations (50) and (51) define a set of plane-waves 
impinging on/emerging from the crossing point r = 
(— £f£b + £o,£f£b — Co) 01 the unperturbed chiral edge 




FIG. 3. (Color online) 2D color plots of |^^(t-)| 2 as ob- 
tained from Eqs. (461 and (471. We remind the reader that 



these solutions refer to the case of linear confinement poten- 
tials, which obey the smooth-confinement condition ( 23 1 . The 



horizontal and vertical axes denote Sx/£b and 5y/iB, respec- 
tively. The numerical results in this plot refer to the case 
in which an electron enters the system from the west in the 



B layer- 



ed = and c 



-nV/2-<£f 



/ 7 r(; 7 2 



The 



left plots refer to the B-layer wavefunction tp^ (r), while the 



(48) right plots to the T-layer wavefunction ip^ (r). The upper 



panels refer to 7 = —0.4, while the lower panels to 7 = — 1.4. 
We have taken £0 = in Eqs. (46 1 and (47 1, consistently with 



the smooth-confinement condition ( 23 1 



modes — the logarithmic phase terms being irrelevant cor- 
rections when compared to the longitudinal phase depen- 
dence of 4> B fe ( ) (r) and (f> T Q fc(0 ) (r). For instance, set- 
ting d = and c = — e^^^ 2 / 2 ~^ /jT(ij 2 ) we obtain 



^ (r) ~ [t Q(Sx) + e(-5x)} B Q fc(0) (r) , 



ip^oir) ~ r 9(Sy) (t> TOJc (o)(r) 



(52) 



which describe a scattering event where an incoming 
wave from the west hand side of the B layer — represented 
by the component ^3 (r) = 0(— 8x)<f> B Q fc (o)(r) — splits 

into a transmitted wave ip^ ut \r) = 6(fe)0 B fe (o) (r) 
propagating on the east hand side of the same layer, and 
into a deflected wave tp^ ut \r) = Q(Sy)4> T . <o) (r) prop- 
agating along the south-north direction in the T layer — 
see Fig. [2] The corresponding transmissivity t and reflec- 
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FIG. 4. The transmission \t\ 2 and reflection |r| 2 probabilities 
of the linear confinement problem as obtained from Eqs. ( 53 1- 



( 54 1 are plotted as functions of I7I . The dashed line represents 

Fr + = i. 



tivity r are determined by the parameters 



7 r(i 7 2 ) 



(53) 
(54) 



which fulfill the normalization condition |r| 2 + \t\ 2 
1 thanks to the identity |r(i7)| 2 = 2tt (e^ 2 



T77 2 



see 



Fig.g 



Notice, in particular, that in the weak tunneling 7 — > 
limit, the reflectivity r in Eq. (54) reduces to the value 



V 27T7 = i e 



71" 1_ bA S as 



(55) 



in perfect agreement with the result based on the Born 
approximation. The latter is obtained from Eq. (21) by 
taking <fi B fc < ) (r) and </> T fc (o> (r) as input and output 
states. 

In a similar fashion, setting b = and a = 

e -^7 2 /2/ 7 r(_j 7 2) ; Eqs _ and O) can be used t0 



describe the scattering event "complementary" to ( 52 ) 
where the incoming wave V't ( r ) = 0(— Sy)(p T k w (r) 
is reaches ro from the south-north direction in the T layer 
and gets partially deflected in the B layer. In this case 
we get 

^ Q (r)~-r* Q(Sx) BAfc <o,(r), 

4«> (r) =s [i*0(<yy) + e(-* V )] T O fe (o> (r) , (56) 

which, together with Eq. (56), ensure the unitarity of the 
mapping [ip^ n] '(r) ,4™' V)] ^ [V>b i W'V't '(»*)]• 



IV. EXACT SOLUTION FOR GENERAL 
POTENTIALS IN THE SMOOTH 
CONFINEMENT LIMIT 

In this Section we address the case of confinement po- 
tentials which are smooth, i, e. they satisfy Eq. ( 23 ) , but 
not necessarily linearly-dependent on the transverse co- 
ordinate u a , as in the previous Section. In this case we 
need to solve Eq. (29) with e CTn fe as in Eq. (25). 



The idea we are going to exploit is to focus the atten- 
tion on an interval of values of k around which the disper- 
sion relation can be linearized. Specifically, we assume 
that around point A in Fig. [l] the following expansion 
holds 



kg +(k- k F )hv F + 



(57) 



for k e (k F - Afc/2, k F + Afc/2). k F is the Fermi mo- 
mentum of the system associated with the energy e, i.e. 
£<r,n,k F — e, while v F is the group velocity of the modes, 

~dk 



dV a (kil) 



k—k^ 



dk 



(58) 



k—kp 



which, for the sake of simplicity, we assume to be layer 
independent. Notice, however, that, in principle, for each 
n one has specific values of k F , Ak, and v F . 

For this purpose, in analogy with the layer field oper- 
ator Vt^r) in Eq. (JsJ> , we now define the layer-resolved 
edge field 



dk 



& 9k Ccr,n,k : 



with re = x (rx = y) for the B (T) layer and 

(A - fc F ) 2 " 



g k = exp 



2Afc 2 



(59) 



(60) 



The introduction of the kernel gk has two consequences: 
i) it forces k to lie in the interval where (57) holds; it) it 



induces a coarse graining in the position resolution of the 
fields, which now obey fermionic anti-commutation rules 
only beyond a length scale 2ir/Ak, 

x S Ak (r a - r' ), (61) 

with 5 Ak (x) = Afc(2V7i : )" 1 exp(-a; 2 Afc 2 /4). 

The edge field ia,n{r) in Eq. (59) can be equivalently 
expanded by using the field 7(e), which satisfies the eigen- 
mode equation (15), 



(62) 



By exploiting Eq. (18), which connects 7(e) to the un- 
perturbed fields c a n kl and which can be inverted as 



E 



* (e) 

a,n,k 



7(e) 



(63) 



we can identify, in complete analogy with the layer wave- 
functions Eq. (19), the coarse grained edge wavefunc- 
tions, 



dk 



9k a k 



(64) 



In the smooth-confinement approximation, Eq. ( 23 1 , 



the complex amplitudes o$ k satisfy Eq. (|29|) . The in- 



troduction of the distribution effectively linearizes the 
dispersion e atn> u for those k which belong to the interval 
(k-p — Afe/2, fcp + Ak/2). It then follows that the wave- 
functions Casifrv) satisfy the following equation 



'7 



'2tt 



dy £Uy)K(x,y), 



(65) 



{k F + id y )^ n (y) 



?"7 



dx&Jx)K*{-y,-x), 



where the kernel K(x, y) is obtained under the condition 
AMb 3> 1 and reads 



K(x,y) 



1 



exp 



(y 



Ml) 2 



2Ak 2 ei 



.xy 



(66) 



By further tightening the condition Ak£s 3> 1 we can 
cast the system of equation for the wavefunctions £^n(r<j) 
in the form 



(A*- 
(fc F 



id x )$ n (x) 
id y )£Uy) 



{-i) n i^U-y/il) 



(67) 
(68) 



The solutions of these equations are given by 
Eqs. ( A7|A10 ), which, by neglecting the logarithmic 
phase exp(±«7 2 In |r CT /^B|) with respect to the linear in- 
crease, read 



(e) (y) 



e ifex e i(W B ) 2 [ a 9^) + b Q(- x )] , (69) 
e^y[ce(y)+dO(-y)]. (70) 



The coarse-graining procedure has averaged out all the 
details below a scale 2ir/ Ak, which is much smaller than 
the magnetic length ^, which in turn is the smallest 
length scale in the problem. Eqs. ( 69|70 ) are in com- 
plete analogy with Eqs. (50 51 1. We can then extend 



the results for the reflection and transmission amplitudes 
Eqs. (53 54 1 derived in the special case of linear confine- 



ment to the general case of smooth confinement. 



in the bottom layer, and "TN", which propagates from 
south to north in the top layer. These two modes cross 
at position A in Fig. [T] In particular, we calculated the 
corresponding transmission amplitude t, Eq. (53), and 
the reflection amplitude r, Eq. (54). 



The scattering at the other crossing points B, C, and 
D can be characterized in a completely analogue way and 
we can write the scattering matrix Sj characterizing each 
beam-splitter i = A, B, C, and D as 



(71) 



with ti and Ti given by t and r, Eq. (53 54), respectively, 
for nominally equal beam-splitters. 

The full scattering matrix S that characterizes the in- 
terferometer response can now be obtained through a 
concatenation procedure^, staring from the scattering 
matrix Sj of the elementary beam-splitter. For instance, 
the amplitude S21 for scattering from contact 1 to contact 
2 — see Fig. [l] — can be obtained as a geometric series that 
sums all the possible paths an electron can take before 
exiting from 2, and for nominally equal beam-splitters 
takes the form 



S21 = r - e iVd ^ 



\t\ 2 (r*) 



* \3 



1 — (r*) 4 e i¥,d >- n 



(72) 



where the phase tfdyn = 4fcpW is the dynamical phase ac- 
quired around the perimeter of the square defined by the 
corners i = A,B,C, and D in Fig. [l] The transmission 
probability T 2 i = IS21I 2 can be written as 



T 2X =R 



1 



R(l-R) 



1 - R 2 e^ 



(73) 



where R = \r\ 2 and ip — fdyn + 4(fcp^s) 2 , up to an offset 
given by the argument of T(—ij 2 ). We notice that the 
phase 4(fcF^s) 2 is equal to the Aharonov-Bohm phase 
27t$/$ ($q = hc/e being the flux quantum), where $ is 
the flux through the area of the square ABCD defined by 
the edge state crossing points of coordinates 2Xj = i/cp^fg 
and yi = ikpig. 

T21 is zero for R — (no tunneling) and one for R = 1, 
(all the current injected in contact 1 is totally drained 
at the first beam-splitter at point A). Analogously, one 
can obtain all the other transmission and reflection am- 
plitudes. 



VI. DISCUSSION AND CONCLUSIONS 



V. INTERFEROMETER SCATTERING 
MATRIX 

The full response of the interfero mete r can be studied 
in the Landauer-Biittiker formalis m 1 32 * 33 ! by means of the 
scattering matrix S. In the previous Sections we have fo- 
cused the attention on the mixing between two particular 
edge states, "BE", which propagates from west to east 



We have proposed a strategy for the implementation 
of electron interferometry which is based on quantum 
Hall bilayers. Our approach exploits uniform inter-layer 
tunneling to mix edge modes. While tunneling between 
the two layers is present in the whole region where top 
and bottom Hall bars overlap (see Fig. [I]), the quasi one- 
dimensional character of edge states implies that inter- 
edge tunneling is active only in a region of linear size 
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of order of the magnetic length around the geometrical 
crossing points of edge states belonging to different lay- 
ers. 

We have demonstrated that the scattering problem as- 
sociated with the quantum Hall bilayer depicted in Fig. [T] 
can be solved exactly provided that the confining poten- 
tials defining the two Hall bars obey the "smooth confine- 
ment" condition (23). The analysis of this Article relies 
on a single-particle picture and therefore neglects inter- 
layer electron-electron interactions. Fascinating many- 
body phenomena including spontaneous inter-layer co- 
herence and Coulomb drag effects are expected to occur 
in QH bilayers^. These effects are particularly spec- 
tacular in the regime in which A§as is a negligible en- 
ergy scale. In this Article we have analyzed the oppo- 
site regime. Electron-electron interactions in the regime 
Asas C e 2 / ((d), where e is a material parameter (e ~ 13 
for GaAs), might have a non-trivial impact on the be- 
havior of the proposed beam splitter. We hope to tackle 
this intriguing regime in a forthcoming publication. 

Turning to interferometric implementations, the key 
advantage of the proposed setup is that it bypasses 
the use of quantum point contacts to mix edge modes. 
"Non-simply connected circuits" are not necessary for the 
achievement of our interferometer. This is expected to 
allow the design of smaller and more advanced coherent 
circuits in the quantum Hall regime. 

Our quantum Hall bilayer beam splitter is expected to 
display a behavior that closely mimics the one of a semi- 
transparent mirror in conventional optics. One drawback 
of our proposal is linked to the fact that in bilayers real- 
ized in semiconductor hetero junctions the magnitude of 
the interlayer tunneling Asas is mostly determined by 
the heterostucture design (and thus fixed by the growth 
procedure) . On the other hand, AgAS can be tuned by an 
in-plane magnetic fielcpSI, which, for sufficiently thin sam- 
ples, does not introduce additional orbital effects. This 
implies a certain degree of tunability of the beam split- 
ter transmittivity. Another degree of tunability, which 
we have not exploited in the present Article, is offered by 
the relative angle 9 between edge modes living in different 
layers (in this Article we have analyzed only 9 — tt/2). 

A very interesting perspective is also constituted by the 
possibility to implement a beam splitter in the fractional 
quantum Hall regime, which has been historically hard. 
While it is not obvious at all how to generalize the current 
theory to such regime, a similar phenomenology might be 
expected and might shed new light on transport phenom- 
ena involving fractional edge modes. Last but not least, 
we would like to mention that coherent beam splitters 
may also be realized in vertical heterostructurea 3 -^ com- 
pri sing g raphene as well as other two-dimensional crys- 
t ai d 36 ' 37 1 such as MoS 2 , h-BN, etc. The theoretical com- 
plication posed by these structures, however, is highly 
non-trivial since a tunneling Hamiltonian as simple as 
the one in Eq. ^ does not apply. Inter-layer transport 
in these systems is very interesting^ and currently far 
from being completely understood. 
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Appendix A: Derivation of the solutions (41 1 
and d42l) 



A convenient way to write Eq. (38) is by introduc 



ing the dimensionless variable k = Mb and the complex 
functions /b(«0 and /t(k) defined implicitly by the ex- 
pressions 



» - 



fc = h{M B ) h 



» - 



fc = {-i) n h(kl B ) l B ■ 



(Al) 
(A2) 



Via Fourier transform Eq. ( 38 ) can then be expressed as 
( K F + i£) /b(«0 =7/t(-k) , 

(A3) 

(k f + k) /t(-k) = 7 /b(«) , 
or, equivalently, as 

(k f - «) /b(«) = 7 /t(«) , 

(A4) 

(kf + i£) 7 

where up stands for the dimensionless parameter k f = 

kp" £bi while we used the symbol /(k) to represent the 
Fourier transform of the function /(«;), i.e. 

1 

/2tt 



/(*) 



(A5) 



By a close inspection of these equations, one notices that 
they fulfill certain symmetries. In particular they are 
invariant in form by replacing /b(k) with the complex 
conjugate of /t(k) and vice- versa, i.e. 



(/bW,/tW)^(/tW-/bW) 



(A6) 



Furthermore it is easy to check that the solutions /b(^), 
/t(k) for the case ftp 7^ can be written as 



; B ( K ) = e i ^+^;£V + KF ) 

f T ( K ) = e ilt ^-^f^ ) (K-K F ) 
or, equivalently, 



"""72V -«f), 



/b(«) = e 

; T ( K ) = e^ K ;£V- KF ), 



(A7) 
(A8) 

(A9) 
(A10) 



with /£V),/tV) solving Eq. Q [or |A4| ] for k f = 
0. Explicit expressions for the latter can then be obtained 
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by observing that (k) fulfills the following (complex) 
Euler differential equation, 



-7 2 i| 0) (-) 



(All) 



[this can be verified via elementary algebraic manipula- 
tion of Eq. (A3)]. Upon integration, we thus obtain 



/< 0) (k) = [a 0(«) + b ©(-«)] | K |-^ 2 , (A12) 



and, via the second line of Eq. (|A3|) , 

»(«) 



4 0) ' 



[a9(-K) 



6 0(«)] |k| 



-*7 



(A13) 



which, together with ( A8 ) , gives us ( 42 1 when replaced 
into (A2|. In these expressions Q(«) is the Heaviside 



step function, while a, 6 are complex coefficients that 
are fixed by the boundary conditions fulfilled by (k) 



and ftp (k). One might notice that the solutions (A12) 



and (A13) are well defined everywhere but for k = 
where they are unstable due to the fast oscillating term 
|k| -17 " and, in the case of fj?'(n), due to the presence 
of the factor 1/k [notice however that the presence of 
"stronger" discontinuities (e.g. Dirac delta or differen- 
tial Dirac delta contribution terms) can be excluded by a 



close inspection of the differential equation (All |]. These 



irregularities cannot be avoided and need to be properly 
accounted for when producing the functions (n) and 

/g (k) via Fourier (or inverse Fourier) transformation. 
The correct prescription is obtained by adopting the fol- 
lowing integral expressions 



dn' 



r( ?7 2 



-1 ^ihtK 



0(«) 



lim 



r( i7 2 + 5) 



dn' 



s^o+ (-iK) in / 2 + 5 
^ 2/2 + Q(-k) e^ 2/2 

r(; 7 2 + s) 



lim 

<5->i- 



r(i 7 2 + 1) 

| K |i7 2 +l 



(-in) i 'r 2 + s 
0(k) e~^ 2 /2 _ (. 



(A14) 



(A15) 



^7/2 



which hold for all k ^ from the identity 3.381(7) 
of Ref. 39 via analytic continuation [here T(- • • ) is the 
Euler function]. 



Accordingly one can easily verify that the following 
identity holds 



^ 0) (-) 

/4 0) («) 



[c Q(k) + d 0(-k)] , (A16) 

[c6(k) +d0(-/s)] |«|^ 2 , (A17) 



with c and d begin complex parameters which depend 
from a and & via the linear constraints ( 45 1 . To check 



the consistency of the procedure observe that following 
the same simple manipulations that led us to Eq. (All I, 
from Eq. (|A4|) we get 



d a 



(0) 



K dn /t 



( K )= J7 2 /4 0) ( K ) 



(A18) 



which indeed admits the functional form of Eq. ( A16) as 
most general solution. 
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